Anomalous transport in biased ac-driven Josephson junctions: Negative conductances 
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We investigate classical anomalous electrical transport in a driven, resistively and capacitively 
shunted Josephson junction device. Novel transport phenomena are identified in chaotic regimes 
when the junction is subjected to both, a time periodic (ac) and a constant, biasing (dc) current. 
The dependence of the voltage across the junction on the dc-current exhibits a rich diversity of 
anomalous transport characteristics: In particular, depending on the chosen parameter regime we 
can identify so termed absolute negative conductance around zero dc-bias, the occurrence of negative 
differential conductance and, after crossing a zero conductance, the emergence of a negative nonlinear 
conductance in the non-equilibrium response regime remote from zero dc-bias. 

PACS numbers: 74.25.Fy 85.25.Cp 73.50.Td 05.45.-a, 



I. INTRODUCTION 

The Josephson junction system constitutes a beautiful 
paradigm of a nonlinear system exhibiting most inter- 
esting classical and quantum phenomenal This system 
in addition offers a rich spectrum of beneficial applica- 
tions. For example, a prominent application of a Joseph- 
son junction set-up relates to the definition of the volt- 
age standard. Moreover, practical devices based upon 
the characteristics of a Josephson junction are instru- 
mental for high speed circuits: They can be designed to 
switch voltage within a few picoseconds. Their typical 
attribute of a low power dissipation proves them service- 
able in high-density computer circuits where the resistive 
heating limits the applicability of conventional switches. 
More recent applications refer to quantum computing de- 
vices where Josephson junction set-ups can store single 
units of information (qubits)^ or with appropriately en- 
gineered coupling among various units, they serve as an 
architecture for the processing of quantum information. 3 
Yet, there remain still new phenomena to be uncovered 
for this system which in turn carry the potential for new 
applications. Of special interest are novel transport phe- 
nomena in presence of external ac- and dc-forcing. In this 
context, the phenomenon of absolute negative conduc- 
tance (ANC) plays a particularly intriguing role. Here, 
ANC means that the system's response is opposite to a 
small external bias. In the present context, a small, say, 
positive dc-current generates a negative voltage across 
the junction. This constitutes no contradiction to ther- 
modynamic laws because it occurs in presence of simulta- 
neous acting non-equilibrium perturbation, the ac-drive. 
This phenomenon is known to emerge within a quantum 
mechanical setting in presence of tunneling processes; re- 
markably, however, it also has been demonstrated on a 
classical level in stylized, spatially extended ratchet-like 
systems without reflection symmetry4^i&£ 

The phenomenon of ANC has already been experi- 
mentally observed in p-modulation-doped GaAs quan- 
tum wells 8 and also in semiconductor superlattices, oc- 
curring therein as a genuine quantum phenomenon^ In- 



terestingly enough, this very ANC-phenomenon has been 
reported first in a recent work by us in Ref. J^Oj for an 
ac- and dc-driven Josephson junction. Notably, however, 
the study in Ref. |RJ identifies this phenomenon within 
its classical operation regime for an inherent reflection- 
symmetric system. In extending our previous study in 
Ref. [1(1 we (i) provide further details, (ii) explore even 
wider regimes in parameter space and (iii) identify ad- 
ditional novel response regimes. In contrast to a related 
studyjii where ANC is investigated in an underdamped, 
deterministic chaos regime, here we emphasize the role 
of thermal noise and noise-induced nonlinear response 
phenomena. Moreover, we put forward a study of an 
accessible, optimal parameter regime towards the objec- 
tive for an experimental verification of our findings. The 
underlying dynamics of this driven Josephson junction 
can conveniently be described by the model of a resis- 
tively and capacitively shunted junction in terms of the 
so called Stewart-McCumber modeli^ 

The layout of the present work is as follows: In Sec. 
II, we present the classical Stewart-McCumber model of 
the Josephson junction. Next, in Sec. Ill, we elucidate 
various regimes of anomalous transport behavior, such 
as ANC, negative differential conductance (NDC) and so 
termed negative- valued nonlinear conductance (NNC). 

In Sec. IV, we work out the optimal regimes for the 
phenomenon of the negative conductance. In Sec. V, we 
elaborate on the experimental feasibility for an in situ 
confirmation of our diversified novel theoretical predic- 
tions. Sec. VI provides a summary and some conclu- 
sions. 



II. THE STEWART-MCCUMBER MODEL 

This model describes the (semi)-classical regime of the 
voltage-current characteristics of a Josephson junction. 
The model involves the Josephson supercurrent charac- 
terized by the critical current Iq, a normal (Ohmic) cur- 
rent characterized by the normal state resistance R and a 
displacement current accompanied with the junction ca- 
pacitance C . The ubiquitous thermal equilibrium noise 
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consists of Johnson noise associated with the resistance 
R. The dynamics of the phase difference 4> — <f>{t) across 
the junction is then described by the following nonlinear 
equation of motion, e.g. see in the review Il2l. 
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Herein, a dot denotes the differentiation with respect 
to time t, Id and I a are the amplitudes of the applied 
dc- and ac-currents, respectively, Q is the angular fre- 
quency of the ac-driving. The parameter k B denotes the 
Boltzmann constant and T is the temperature of the sys- 
tem. Thermal equilibrium fluctuations are modeled by 
J-correlated Gaussian white noise £ (t) of zero mean and 
unit intensity, i.e., < £(t)£(s) >= S(t — s). 

The limitations of the Stewart-McCumber model and 
its range of validity are discussed e.g. in Sec. 2.5 and 
2.6 of the comprehensive review paper by Kautz^ In 
particular, we thus work within the small junction area 
limit and in a regime where photon- assisted tunneling 
phenomena do not contribute. Throughout the following 
we shall adopt the dimensionless form of Eq. (TTJ) from 
Refs. [l| and[T3, namely, 
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+ cr-^j + sin(0) = i + ii cos(fiit') + V2aD T(t'), 
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where the dimensionless time t' = f/ro and the charac- 
teristic time To = 1/ujp. The Josephson plasma frequency 
uj p = (l/h)\/8EjEc is expressed by the Josephson cou- 
pling energy Ej = (h/2e)Io and the charging energy 
E c = e 2 /2C. The 'friction' coefficient a = t /RC is 
given by the ratio of two characteristic times: To and the 
relaxation time r r = RC. This dimensionless friction 
parameter a measures the strength of dissipation. The 
amplitude and the angular frequency of the ac-current 
are i\ = I a /Io and fix = flro = fl/to p , respectively. 
The re-scaled dc-current reads io = Id/Io, the re-scaled 
zero-mean Gaussian white noise T(t') possesses the auto- 
correlation function (T(t')T(u)) = 5(t' — u), and the noise 
intensity D — k B T/Ej is given as the ratio of two en- 
ergies, the thermal energy and the Josephson coupling 
energy (corresponding to the barrier height). A different 
scaling procedure, being more familiar within the Brow- 
nian motor community, is detailed in Appendix. 

The most important characteristic for the above sys- 
tem is the current-voltage curve. To obtain it, we numer- 
ically integrated Eq. ([2]). For long time (to avoid initial 
conditions and transient effects), we calculated the sta- 
tionary dimensionless voltage 



'dt' 



(3) 



over one cycle period of the external ac-driving. 
stationary physical voltage is then expressed as 

V = — - v. 

2e 
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where the brackets denote an average over both, all the 
realizations of the thermal noise and a temporal average 



Strictly speaking, the zero temperature limit, D = 0, 
should not be considered within the framework of the 
Stewart-McCumber model; because at zero temperature 
the quantum dynamics comes into play. However, in or- 
der to explain peculiar properties of the system at non- 
zero temperature and/or to seek optimal regimes, it is 
insightful nevertheless to consider the deterministic dy- 
namics, i.e. the zero noise case, D = 0. Then, the role of 
initial conditions are also relevant. This is particularly 
so if several attractors coexist. In such a case, an addi- 
tional average over initial conditions must be performed: 
We have chosen initial phases 4>q that are equally dis- 
tributed over one period [0, 2ir] and (dimensionless) ini- 
tial voltages vo equally distributed in the range [—2,2]. 
Further details of how to treat the deterministic case are 
presented e.g. in Ref. [l3|. In the case of several attrac- 
tors, it means that the average is over attractors, whose 
weights are proportional to corresponding basins of at- 
tractions. However, it is worth to stress that the case 
D = generally is not equivalent to the limit D — > 
performed in the non-zero temperature, Fokker-Planck 
case, see the studies in Ref. [l^. 

Physical systems described by Eq. |J5J are widespread 
and well known. An example is a Brownian particle mov- 
ing in the spatially periodic potential U{4>) — U(<fi + L) = 
— cos((/>) of period L — 2tt, driven by the time-periodic 
force and a constant forced Then, the variable (f> corre- 
sponds to the space coordinate of the Brownian particle 
and the ac- and dc-currents plays the role of driving forces 
on the particle. Other specific systems include rotating 
dipoles in external fields , 16 ' 17 superionic conductors^ or 
charge density waves^ to name just a few. 



III. ANOMALOUS TRANSPORT BEHAVIOR 

The deterministic dynamics corresponding to @ en- 
compasses a three-dimensional phase space, namely 
{4>,d(f)/dt' and contains four parameters, reading 
{ct, «o, «i, Therefore, its dynamics is able to exhibit 
an extremely rich behavior in phase space as a function 
of the chosen system parameters! 12 i 20 For example, its 
dynamics features harmonic, subharmonic, quasiperiodic 
and also chaotic types of behavior; for further details we 
refer the readers to the review in Ref. fl2l . At non-zero 
temperature, D > 0, thermal fluctuations lead to diffu- 
sive dynamics for which stochastic escape events among 
possibly coexisting attractors are possible. 21 Moreover, 
when thermal noise is acting, the system dynamics is con- 
stantly excited away from stable trajectories; it thus can 
explore the whole phase space. In some cases the prob- 
ability distribution can concentrate on regions in phase 
space which do not coincide with stable trajectories. The 
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prominent example of such phenomena is a classical ex- 
citable stable point, which under influence of noise can 
produce a "limit cycle"-like probability distributions. 22 
This in turn can imply drastic consequences for the trans- 
port properties. 

From the reflection symmetry <p — > —0 of the poten- 
tial U(4>) — — cos((/>) and the time-reflection t 1 — > — t' 
of the ac-driving, it follows that the average voltage is 
strictly zero if the dc-current is zero (io = 0). If a nonva- 
nishing dc-current is applied, however, (i.e. io ^ 0), the 
above mentioned reflection symmetries are broken and a 
non-zero voltage typically emerges. Since the dynamics 
determined by Eq. ([2]) is nonlinear and the system is 
multidimensional, it should not come as a surprise that 
the current-voltage characteristics is typically nonlinear 
and often depicts a non-monotonic function of the sys- 
tem parameters. Nevertheless, some most non-intuitive 
behaviors still remain to be unraveled in parameter space 
which seemingly have escaped previous detailed investi- 
gations for this archetype system 1 ^. 
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FIG. 1: The diagram depicts the dependence of the volt- 
age versus the variation of the driving amplitude i\ of the 
ac-current for a positive dc-current set at io = 0.0159. Re- 
maining parameters are chosen as a — 0.143 , fii = 0.78 and 
temperature D — (dotted line, blue online, left scale) and 
D = 10 -3 (solid line, black online, right scale). 



The current-voltage curve v = w(«o) is a nonlinear func- 
tion of the dimensionless dc-current strength io- Upon 
inspecting the symmetries in the equation of motion 
we find that this function is odd in the dc-bias io, i.e. 
v(— io) — —v(io). Typically, the voltage is an increasing 
function of the dc-current. Such regimes correspond in 
parameter space to a normal, Ohmic-like transport be- 
havior. More interesting are, however, those regimes of 
anomalous transport, exhibiting (i) an (absolute) nega- 
tive conductance near zero dc-bias, (ii) the regime of a 
negative differential conductance, and (iii) after cross- 
ing zero conductance, a negative-valued conductance in 
a nonlinear regime displaced from zero dc-bias io = 0. 

The numerical analysis of the system dynamics de- 
picts that anomalous transport occurs in a parameter 
range where the driven system dynamics is strongly non- 



linear. Although there is no obvious direct connection 
to chaotic properties of the system dynamics we have 
found that regimes of anomalous transport typically ne- 
cessitate also a chaotic dynamics. In particular, for the 
regime a = 0.143 and fti = 0.78 presented in Fig. Q] the 
numerical study shows that the first bifurcation cascade 
at D — leading to chaos occurs near i\ = 0.477, c.f. the 
dotted line in Fig. [TJ while the regime of an anomalous 
transport behavior starts out near i\ ~ 0.637. Not un- 
expectedly, no anomalous transport occurs in the regime 
of approximate linear dynamics. 



A. Absolute negative conductance 

In transport theory, the Green-Kubo linear response 
regime plays an important role. It allows one to obtain 
linear transport coefficients. For the system described by 
Eq. (0), there are regimes, where for sufficiently small 
values of the dc-bias «o, linear response theory holds. 
This is characterized by the relation 23,24 



(5) 



It defines the linear transport coefficient Vl = 
r L(c, ii,fl\,D) which does not depend on io and is called 
the static resistance or = 1/tl, is the conductance in 
the linear response regime. Since the full response ex- 
hibits also a non-linear behavior, it is necessary to recall 
that the resistance is defined in the limit of a small 
dc-current, r l = limj — >o[ v (io)/io]- The case > cor- 
responds to the normal transport behavior or the Ohmic- 
like regime. The case tl < amounts to the situation 
when the voltage assumes the opposite sign of the dc- 
bias. It is termed here absolute negative conductance 
(in analogy to the absolute negative mobility of driven 
Brownian particles^&i) . In Fig. [5] we exemplify this 
situation. Indeed, for dc-current values varying between 
io € (—0.04, 0.04), the voltage assumes a sign opposite to 
io- The driving amplitude of the ac-current at i\ = 0.73 
corresponds to the second negative-valued minimum of 
the averaged voltage shown in Fig. [1] for a temperature 
D = 10~ 3 . The first negative- valued minimum of v ver- 
sus i\ in Fig. Q]has been explored in greater detail in our 
previous work in Ref. [ToL Let us emphasize here that, 
in distinct difference to this case with non-zero thermal 
noise D ^ 0, in the limit of vanishing noise (i.e. when 
D = 0) the averaged voltage is in fact identically zero at 
the ac-driving strength i\ = 0.73. This implies that this 
counter-intuitive phenomenon of ANC is solely induced 
by thermal equilibrium fluctuations. For the record, there 
exist also linear response regimes where even in the de- 
terministic case, the voltage can become negative for a 
positive dc-current (and vice verse) i 10 ' 1:L From Fig. [TJ one 
can observe that there exist several parameter windows 
depicting such a noise-induced ANC-behavior: the con- 
ductance <7l = 5l(c, ii, &i,D) exhibits sign changes mul- 
tiple times upon increasing the ac- amplitude strength i\. 
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Between such ANC windows with g^ < 0, normal trans- 
port regimes with > occur. Moreover, we tested that 
the ANC windows remain stable upon a small variation 
of the remaining parameters. 



B. Negative differential conductance 

If the voltage v is not a monotonic function of the 
dc-current io, the differential conductance can assume 
negative values. The differential (or dynamic) resistance, 
defined by the relation^ 4 - 



Td (io) 



dv{ig) 
din 



(6) 



may therefore assume negative values within some inter- 
val of dc-bias values in. In Fig. [21 we depict two examples 
of such behavior: the negative differential conductance 
(NDC) go{io) — l/ro(io) is detected for the dc-current 
i e (0.0796,0.086) and i e (0.132,0.166), note the two 
inset panels in Fig. [21 To the best of our knowledge, 
this typical NDC voltage-current characteristic shown in 
Fig. [3 has not yet been experimentally reported for sym- 
metrical systems, such as the system in Eq. @. For 
asymmetric ratchet systems, this effect was described in 
Ref. [H 

We note that both effects of absolute and differen- 
tial negative conductance are found in the part of the 
£7i — a plane where Oi < 1 and £7i > a. In this re- 
gion Shapiro steps do not exist, i.e. there are no phase 
locked regions which interrupt an otherwise linear current 
current-voltage characteristic. For a detailed analysis of 
the different regions of the Sl\ — a plane with respect 
to the occurence of Shapiro steps see the discussion in 
chapter 5.1 of the review by Kauta 1 -^. In Fig. [5] current 
voltage characteristics for two parameters sets Oi = 0.78, 
a = 0.143 and fii = 7.8, er = 0.143 are displayed. The 
characteristics for the smaller frequency contain several 
bias current intervals with negative conductance, several 
locking regimes but no Shapiro steps. The second cur- 
rent voltage characteristic does have pronounced Shapiro 
steps but no regions with negative conductance. 



tends to zero. In the nonlinear response regime, the volt- 
age can tend to zero even if the dc-current assumes a 
non-zero, finite value. For example, there can emerge 
two positive values of the dc-current {ioi,io2} such that 
in-between the voltage is negative. Then, the nonlinear 
conductance coefficient gN(io) is negative- valued in this 
very interval (ioi, io2); meaning that the voltage takes on 
the opposite sign of the dc-current for a non-zero valued 
dc-bias strength io- Indeed, this situation is presented 
with Fig. [3] The amplitude of the ac-current is set at 
i\ = 0.668; it corresponds to the first negative-valued 
minimum of the averaged voltage depicted in Fig. [TJ For 
small dc-currents, ANC exists (see inset of Fig. [3]) while 
for larger dc-currents io, negative- valued nonlinear con- 
ductance (NNC) results. NNC is a predominantly deter- 
ministic phenomenon which survives in presence of small 
thermal noise. Indeed, we observe in Fig. [3] that for the 
case D = 0, this effect is most pronounced and takes 
place in a wide parameter interval (ioi,io2)- The long- 
time trajectories of the dynamical system ([2]) correspond 
to running states into the negative direction of </>, which 
are of period two and therefore the deterministic aver- 
age voltage is large, cf. Fig. [31 When the temperature 
increases, the interval (ioi,io2) shrinks and the ampli- 
tude of the voltage decreases. Above some temperature 
(D as 8 • 10~ 4 in Fig. [3]), the NNC effect disappears and 
only ANC survives. 

We end this section by a statement that the occur- 
rence of anomalous transport may be governed by dif- 
ferent mechanisms. In some regimes it is solely induced 
by thermal equilibrium fluctuations, i.e. the effects are 
absent for vanishing thermal fluctuations D — 0. In 
other regimes, anomalous transport may also occur in the 
noiseless, deterministic system and the effects increas- 
ingly fade out with increasing temperature*-. Both situ- 
ations are rooted in the complex deterministic structure 
of the nonlinear dynamics governed by a variety stable 
and unstable orbits. 



IV. OPTIMAL PARAMETER REGIMES FOR 
THE OCCURRENCE OF NEGATIVE 
CONDUCTANCE 



C. Nonlinear response regime displaced from zero 
bias io — 0: negative-valued conductance 

In the nonlinear response regime, the "nonlinear resis- 
tance" or the static resistance at a fixed bias current is 
defined by the relatio n 23 i 24 

r N (i Q ) = r N (i ;a,ii,ili,D) = . (7) 

It typically depends on io in a nonlinear and non- 
monotonic manner. The nonlinear conductance refers to 
its inverse, i.e., g_/v(io) = l/r/v(io)- I n the linear response 
regime, the voltage tends to zero when the dc-current 



In our recent work in Ref. [lCj and in the previ- 
ous section, we discussed some fixed parameters for 
{a, Oi} for which the negative conductance (ANC or/and 
NNC) does emerge. Are these two values exceptional? 
To answer this question, we have searched the part of 
the three dimensional parameter space that is specified 
by ac-driving strength {ii S (0,6.37)}, angular driv- 
ing frequency {Oi 6 (0, 1.91)}, friction strength {a E 
(0,0.796)} in order to locate regions of negative conduc- 
tances in the deterministic case, when D = 0. We set 
the dc-current to a small value io = 0.0159 and calcu- 
lated the value of the current for many randomly chosen 
parameters it, f2i and cr^l The inspection of the results 
revealed that there are values of friction a for which re- 
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gions of negative conductance are most prominent. One 
such value is a = 0.191 for which the negative conduc- 
tance is most pronounced in a relatively large domain of 
parameter variation with relatively large values of the di- 
mensionless voltage. Therefore we performed a more ac- 
curate search for the section a = 0.191 and analyzed the 
two-dimensional parameter domain {i\ S (0,6.37),ili € 
(0, 1.91)}. The results are depicted with Fig. g) The 
points in Fig. [4] which are indicated by different grey 
scales correspond to values where v(i — 0.0159) < 0: 
the voltage assumes thus the opposite sign of the dc- 
current. Strictly speaking a negative voltage at the fi- 
nite bias io — 0.0159 does not necessarily imply ANC. 
However, in most of cases, under these conditions ANC 
is observed in the presence of small thermal noise. We 
have also found, that the noise induced ANC typically 
occurs also for parameters which differ only slightly from 
those where v(io = 0.0159) < 0. Thus, this systematic 
analysis provides insight into the structure of the param- 
eter space and proves useful for designing corresponding 
Josephson junction experiments. Finally, using exten- 
sively this technique we are rather convinced that no re- 
maining regimes of the negative conductance are likely to 
emerge in this system than those already depicted in Fig. 
IH The grey scales in Fig. |4] represents various regimes 
of the ratio v/Q,\. One can distinguish two features in 
the parameter space: a stripe-like structure for a broad 
range of parameters and one pronounced region, being 
zoomed in panel (b). The stripe- like structure suggests 
that in an experiment, one will observe with large prob- 
ability negative conductance upon variation of a single 
parameter. In this context the region shown in panel (b) 
is especially interesting because the negative conductance 
depicted there exhibits a relatively large robustness with 
respect to a variation of the system parameters. 

In the deterministic case, we identify three classes of 
trajectories which generate negative conductance. This 
is marked by different grey scales in Fig. |4j Dark grey 
corresponds to chaotic trajectories. For this regime of pa- 
rameters the conductance is negative also for small non- 
zero temperatures. Light grey corresponds to periodic 
orbits. For this regime of parameters the conductance 
is negative also for small non-zero temperatures. Black 
regions also correspond to periodic orbits; for this regime 
of parameters the negative conductance however quickly 
diminishes as the temperature is raised. 

From the experimental view-point, the most interest- 
ing regimes are those where the negative conductance is 
most robust within some finite temperature interval. In 
the case presented in Fig. [U this corresponds to the dark 
and light grey regions. Therefore, we set e.g. fli = 1.1, 
cf. panel (b) of Fig. [4j in order to observe the variation of 
the negative conductance when passing through the dark 
and light grey regions. Details are depicted in the bifur- 
cation diagram, see Fig. \S[ In the deterministic case, the 
voltage remains almost constant under small variations 
of the amplitude i\ of the ac-current when the param- 
eters belong to the light grey regions, note the voltage 



dependence on i\ around i\ = 2.64 and i\ = 3 depicted 
in panel (a) of Fig. [5] for D = 0. On the other hand, the 
voltage changes irregularly when i\ is changed smoothly 
within the red regions, note the voltage dependence on 
i\ around i\ — 2.78 depicted in panel (a) of Fig. [5] for 
D = 0. In Fig. O the averaged voltage is also shown for 
a small but non-vanishing temperature. When driving 
within the interval a S (2.55,3.15), the voltage is nega- 
tive for D = 10~ 3 and relatively stable with respect to a 
small variation of parameters. The region around i\ = 3 
lends itself as optimal for an experimental verification of 
our predicted anomalous transport features. 



V. EXPERIMENTALLY ACCESSIBLE 
REGIMES FOR ANOMALOUS TRANSPORT 

The results presented above are given for dimension- 
less variables and dimcnsionless parameter values. In or- 
der to motivate experimentalists to test our predictions 
and findings, it is convenient to transform all quantities 
back to their original dimensional values and dimensional 
parameter strengths. There are three important parame- 
ters which characterize the function of a Josephson junc- 
tion. These are the critical current Jo, the resistance R 
and its capacitance C . Three further parameters char- 
acterize the external driving, namely the strength of the 
dc-current Id, the amplitude of the periodically varying 
ac-current I a and its angular frequency Q. Finally, the 
temperature T must be chosen large enough such that 
the junction operates in the semi-classical regime. These 
physical quantities are related to the corresponding di- 
mensionless quantities by the relations: 



V 



hhJp 

~2e 



£1 — UJp^ll , 



1 

RC 



(8) 



with the frequency scale io v given by the characteristic 
plasma frequency, reading 



2 _ 2el 
° p ~ HC 



(9) 



The driving strengths and the actual physical tempera- 
ture read 



I a = Ioh, h = Jo*o, T = 



2ek B 



D , (10) 



being all scaled by the value of the critical current Jo- 

First, one should fix the operational temperature T of 
the classical experimental regime. The last relation in 
(fT0| then yields the strength of critical current io- This 
in turn determines the amplitude strengths I a and Id- 
Fixing the frequency v = Q/2tt of the microwave source, 
then determines via the chosen relevant value of the di- 
mensionless parameter Oj the strength of the plasma fre- 
quency w p . This in turn determines the magnitude of the 
capacitance C via the relation in © and the resistance 
R follows from the last relation in (IS]). Moreover, one 
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should check whether the parameters such chosen obey 
the inequalities: 

hup « Ej, hio p « k B T, (11) 

which taken together guarantee that the junction indeed 
operates in the (semi)-classical regime as presumed with 
the model equation in (p}. These inequalities imply that 
the level spacing of the plasma oscillation is small, both 
compared to the coupling energy Ej and the thermal 
energy ksT. 

To be explicit, we here evaluate some real experi- 
mental circumstances which can present "optimal" con- 
ditions to experimentally verify negative conductance. 
Upon inspection of our dimensionless analysis the follow- 
ing parameter set is an example of the optimal regime: 
i x = 3, fii = 1.11, cr = 0.191 and D = 1CT 3 . For 
a physical temperature of T = 4K, the critical cur- 
rent then is io = 167.8/xA, the amplitude strength be- 
comes I a = 507/iA, the ac-angular frequency emerges as 
fl = 150GHz, Id = 2, 67/iA with the capacitance being 
C = 27.9 pF and the resistance value at R = 1.4S1. Un- 
der these conditions, the absolute value of the voltage 
amounts to V — 3.54^V. 

The NDC regime, which is presented in Fig. [2l is ob- 
served for h = 0.73, Q,i = 0.78, a = 0.143 and D = 10~ 3 
which for a temperature at T = 4K implies a critical 
current Iq = 167.8/iA, I a — 122/iA, ac-angular frequency 
il = 150GHz, and I d varying between € (13.36, 14.43)/^A 
for io G (0.0796,0.086), with the capacitance set at 
C = 13.83 pF and the resistance set at R — 2.63S1. Un- 
der these conditions, the absolute value of the voltage 
approximately reads V = 730/iV. 

VI. CONCLUSIONS 

With this work we took a closer look at the richness of 
anomalous transport behavior occurring in a biased and 
harmonically driven Josephson junction. As it turns out, 
the underlying chaotic dynamics together with the influ- 
ence of thermal noise triggers a whole new variety of un- 
expected transport features. Apart from regions display- 
ing negative differential conductance behavior we could 
identify novel transport characteristics such as noise- 
induced absolute negative conductance near zero bias and 
negative- valued conductance in the strongly nonlinear re- 
sponse regime. Let us summarize these various transport 
phenomena occurring in an ac-driven Josephson junction 
as described by Eq. |T|): 

(1) The dependence of the voltage on the angular driv- 
ing frequency Qi of the ac-current drive depicts 
windows of deterministic and thermally induced 
ANC regimes: ANC appears and disappears as the 
frequency increases (for a preliminary account on 
this effect see also Ref. |28| ). 

(2) The dependence of voltage versus the amplitude i\ 
of the ac-current drive depicts windows of thermal- 
noise induced ANC regimes, cf. Fig. [T] 



(3) The voltage behavior at fixed bias as a function 
of thermal temperatures D exhibits many famil- 
iar features known from the field of Brownian 
motors^ 9 -, such as the occurrence of a voltage re- 
versal versus D, or a typical bell-shaped behavior 
versus noise strength Dj 3 - - 

(4) The voltage as a function of the bias current io 
can exhibit each of the three anomalous trans- 
port features, namely absolute negative conduc- 
tance (ANC), negative- valued nonlinear conduc- 
tance (NNC) and negative differential conductance 
(NDC). 

(5) Reentrant phenomena of negative conductance 
regimes occur as a function of the dc-bias io: Start- 
ing out from zero, the voltage may decrease for in- 
creasing io, reaching a negative- valued local mini- 
mum value which changes over into a local positive- 
valued maximum upon increasing io. Upon further 
increasing io, the voltage starts to decrease again 
into a local, negative-valued minimum, thus ex- 
hibiting NNC. Finally, it increases monotonically 
with increasing io, displaying an almost perfect 
Ohmic-like dependence. 

Our identified novel transport features as presented in 
Fig. Q] and Fig. [5] are also accessible to an experimental 
verification via appropriately designing the experimental 
working parameters for the Josephson system. Here we 
identified such parameter sets in our section V while yet 
a different one has been indicated in our earlier presen- 
tation in Ref. We are confident that our predictions 
will invigorate experimentalists to undertake the experi- 
mental efforts to check our various predictions. 
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Appendix 

In our previous work 10 we used a scaling popular 
among the 'ratchets' community, where the 'coordinate' 
variable x — (f>/2ir is rescaled to the unit interval. Then 
the rescaled form of Eq. (TJJ) reads 

i + 7i + 27rsin(27ra;) = f + acos(us) + ^firfD T(s), (12) 

where the dot denotes differentiation with respect to the 
dimensionless time s = t/r\, where t± — 2tt/lu p . The 
relations between the parameters in two scalings in eq. 
© and |TJ]) are as follows: 

7 = 2na, f = 2nio, a — 2nii, to = 27rf2i (13) 
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Now, the dimensionless velocity v = (dx/ds) and expres- are identical for both scaling procedures, 
sions for the voltage in Eq. (4) and the noise intensity D 
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FIG. 2: The functional dependence of voltage v on the dc cur- 
rent io (i-v characteristic) is compared for two sets of param- 
eters which either display negative conductance or Shapiro 
steps. Panel (a) which corresponds to ac frequency fii = 0.78, 
ac driving strength i\ — 0.73, friction a = 0.143 and noise 
strength D = 10~ 3 displays both absolute negative mobility 
in an interval containing io = and different regions with neg- 
ative differential conductivity, cf. also the respective enlarged 
regions in the two insets. Panel (b) displays the i-v charac- 
teristic for the parameter values of panel (a) (thin solid line) 
and the i-v characteristic for Qi = 7.8, i\ = 73, a — 0.143 and 
D = (thick solid line). Shapiro steps are clearly visible at 
v — kio/a for k = 1/2, 1, 2, 3 (dotted horizontal lines). Away 
from these steps, the i-v charateristic follows the Ohmic law 
v — io/cr (dotted inclined line). Panel (c) presents the i-v 
characteristic for the small driving frequency (f2i = 0.78) for 
positive bias currents io up to io = 0.48 (solid line), the i-v 
characteristic for the same set of junction parameters at zero 
noise (D — 0) (dotted line with wide distance of dots) and 
compares them with a rescaled version v c := 0.1 • v(10 ■ io) 
of the i-v characteristic for the larger frequency such that the 
full range shown in panel (b) is displayed again. Whereas 
the linear overall behavior of the large frequency curve is in- 
terrupted only at a few Shapiro steps, the small frequency 
i-v characteristic is dominated by the nonlinearities also away 
from the locking regions at v — fii/2 and v = 3f2i. 
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FIG. 3: Voltage vs dc-current: The re-entrant effect of the 
negative response (v < 0) is shown at a positive bias (io > 0) 
for the following set of parameters: i\ — 0.668, a = 0.143, 
Qi — 0.78 and four values of temperature D. For small io, 
the absolute negative conductance is observed for non-zero 
temperature and it is solely induced by thermal fluctuations, 
cf. panel (a). For larger io, negative nonlinear conductance 
results as a deterministic effect. The corresponding long-time 
trajectory 4>(t') for io = 0.0796 shown in panel (b) belongs 
to running states into the negative direction of <j> resulting in 
the large average voltage. The value of the amplitude of the 
ac-driving corresponds to the first negative-valued minimum 
of the mean voltage with respect to the ac-current, cf. in Fig. 
□ 




FIG. 4: (color online) The transport properties of the driven 
system in the parameter space {ii , fii } at a representative 
friction value of a = 0.191, dc-bias zo = 0.0159 and zero 
noise strength D = 0. All points in the parameter space, 
where a negative conductance occurs, are marked by symbols 
with different grey scale. The coding corresponds to different 
regimes of values assumed by the ratio v/fli: dark grey (on- 
line: red) denotes the interval (—0.94,0), light grey (online: 
green): ( — 1.26,-0.94), and black (online: blue) corresponds 
to values less than —1.26. Most of the dark grey points cor- 
respond to chaotic trajectories. It turns out that the black 
regimes are most susceptible to thermal fluctuations, meaning 
that the negative-valued conductance rapidly fades away to- 
wards positive values with increasing thermal noise intensity. 
In panel (a) one clearly can distinguish a large compact re- 
gion exhibiting negative-valued conductance, being enlarged 
with the upper panel (b). In regions different from this one 
we find that negative-valued conductance occurs in narrow 
"bands" only. For an experimental realization of negative 
conductance the most promising parameter regimes are those 
regions marked by light grey. This holds in particular for the 
zoomed light grey regimes depicted in panel (b). 




FIG. 5: ( color online) This figure elucidates the transport 
features and the dynamical properties for parameter values 
in space within the section obtained by fixing the 

angular driving frequency at Qi — 1.1, cf. the dotted line in 
panel (b) of Fig. [4] The blue (solid) line in panel (a) depicts 
the dependence of the average voltage on the amplitude of 
the ac-current for zero noise D — 0. The influence of thermal 
fluctuations on the average voltage is shown for two tempera- 
ture values D = 10~ 3 (dashed line) and D = 3 • 1CT 3 (dotted 
line). We do find that a negative- valued conductance indeed 
can survive at small non-zero temperatures. The transport 
properties of the system follow from the underlying complex 
dynamics. In panel (b) we show the bifurcation diagram: the 
Poincare section of the phase in the deterministic system. The 
attractors are color coded according to the sign of the assumed 
corresponding average voltage: blue yields a negative average 
voltage (also pointed by arrows with v < 0), green gives zero 
(also pointed by arrows with v = 0) and red amounts to a 
positive average voltage (also pointed by arrows with v > 0). 



